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Spin relaxation dynamics in rings with Rashba spin-orbit coupling is investigated using spin kinetic 
equation. We find that the spin relaxation in rings occurs toward a persistent spin configuration 
whose final shape depends on the initial spin polarization profile. As an example, it is shown that a 
homogeneous parallel to the ring axis spin polarization transforms into a persistent crown-like spin 
structure. It is demonstrated that the ring geometry introduces a geometrical contribution to the 
spin relaxation rate speeding up the transient dynamics. Moreover, we identify several persistent 
spin configurations as well as calculate the Green function of spin kinetic equation. 

PACS numbers: 72.25.Rb, 71.70.Ej, 71.10.-w 



INTRODUCTION 

The understanding of spin relaxation dynamics in 
semiconductor structures is an active area of research 
related to the field of spintronics [TJ [5] . Previously, the 
spin relaxation dynamics has been investigated mainly in 
infinite two-dimensional (2D) systems with D'yakonov- 
Perel' O H] spin relaxation mechanism (see, e.g., Refs. 
[4HTT|). There are only several examples in the lit- 
erature where the effects of boundary conditions [T^ 
on D 'yakonov-Perel' spin relaxation have been explored. 
These examples include investigations of spin relaxation 
in 2D channels |13l IT¥| . 2D half-space [T5], 2D systems 
with antidots [16], small 2D systems [17], and finite- 
length one-dimensional (ID) structures flF. The main 
conclusion of all these studies is that the introduction 
of boundary conditions results in an increased spin life- 
time. In particular, in Ref. [T7] Lyubinskiy reports on 
a strong suppression of spin relaxation in small 2D sys- 
tems. In Ref. [TB] the present authors demonstrate that 
a persistent spin helix spontaneously emerges in course 
of relaxation of homogeneous spin polarization. As real 
electronic devices are always of a finite size, the under- 
standing of spin relaxation dynamics in reduced geome- 
tries is of a crucial importance. 

In this paper, we present our studies of spin relaxation 
in rings with Bychkov-Rashba |19| spin-orbit coupling. 
The ring geometry is especially interesting since on the 
one hand, the electron space motion in rings is confined 
to a limited space region and on the other hand, the 
ring geometry does not require any boundary conditions 
(such as those derived in Ref. [T2]) for ID transport that 
normally have to be used, for example, to describe spin 
dynamics in finite- length structures [TH]. Starting with 
a spin kinetic equation, we formulate a set of equations 
for spin polarization components on the ring and solve 
these equations in some particular cases. Specifically, 
we find that a homogeneous spin polarization transforms 
into a persistent crown-like spin structure schematically 
shown in Fig. [TJ Moreover, we determine a set of non- 



trivial spin persistent states that are realized at some 
specific sizes of the ring. Finally, we derive expression 
for Green function that can be used to find evolution of 
spin polarization for any initial spin polarization profile. 
All these results indicate an unusual character of spin 
relaxation in finite size structures that is dramatically 
different from the spin relaxation in the bulk. 



SPIN KINETIC EQUATION 

In this section, we introduce the model and derive 
equations governing spin relaxation dynamics in rings. 
The main results of this section are given by Eqs. (|7])-(|9|. 
The derivation of these equations is based on a kinetic ap- 
proach. It is important to mention that Eqs. (|7])-(|9| are 
more general than traditional spin drift-diffusion equa- 
tions [5DH11|. In particular, Eqs. ([7])-([9]) are applicable 
in both diffusive and ballistic regimes at any value of spin 
rotation angle per mean free path. 

Let us consider spin-polarized electrons that can move 
along a ring of a sufficiently large radius r ^ h/p so 




FIG. 1: (Color online) Schematic diagram of spin relaxation in 
a ring. Initially homogeneous spin polarization in z direction 
(along the ring axis) transforms into a persistent crown-like 
spin structure. Note that the spin polar izat ion amplitude 
decreases in this process according to Eq. (21 1. 
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that small-radius corrections |25l E5] to spin-orbit Hamil- 
tonian |jl9j can be neglected. The spin-orbit Hamilto- 
nian is taken in the form Hso ~ a {& x p) ■ e^- Here, 
P — (PxTpy) is the 2D electron momentum operator, m 
is the effective electron's mass, a is the Pauli-matrix vec- 
tor, a is the spin-orbit coupling constant and is a unit 
vector along the ring axis. It can be shown that a mo- 
tion of an electron with a momentum p along a straight 
trajectory is accompanied by a spin rotation with the 
angular velocity ftp = 51p x e^/p, where = 2ap/h. 
The spin precession angle per unit length is given by 
77 = 2amh^^. The ring geometry, however, changes the 
character of spin rotations since the spin rotation axis 
rotates as electron moves along the ring. 

Our consideration of spin dynamics in rings is based on 
the kinetic equation for electron spin polarization (see, 
e.g., Ref. [171 127j ). In quasi-classic approximation such 
an equation can be written as 



d_ 
di 



m 



V 



St{%}, 



(1) 



where Sp(r,t) is the vector of spin polarization of elec- 
trons, and S'ijSp} is the collision integral describing elec- 
tron scattering processes. Eq. (jlj describes change in 
spin polarization of electrons moving with momentum p. 
The RHS of Eq. ([T]) includes two terms describing rota- 
tion of spin polarization with the angular velocity ilp and 
change in spin polarization due to scattering processes. 
In the r-approximation [28] the collision integral is given 

by 



5t{Sp} = — (Sp-(Sp)), 



(2) 



where the angle brackets denote averaging over direction 
of electron momentum. The collision integral ([2| cor- 
responds to the elastic scattering of electrons by strong 
scatterers with a characteristic time t between the col- 
lisions. For ID case the average spin polarization sim- 
plifies to the following expression (Sp) = (S+ -I- S^)/2, 
where S+ and are the spin polarizations of electrons 
moving along the ring in the clockwise (with momentum 
p = mweg), and counterclockwise (p = —mveo) direc- 
tions with the average velocity v connected to the mean- 
free path dhy i — vt. Thus, the kinetic equation ([T|) for 
ID ring of a radius r takes the form of the system of two 
vector equations 



d V d 
dt ^ rde 
d V d 
dt ~ rde 



r!e, xS+-— (S+-S-), (3) 
2t 



-Qe,. xS- - —{S- -S+). (4) 



This system of equations should be complimented by the 
initial conditions for spin polarizations S~^{9,t — 0) and 
S~{9,t = 0) (for clockwise and counterclockwise moving 
electrons) . 



Taking the sum and difference of Eqs. ([3j |4|) we easily 
obtain 
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(5) 
(6) 



where the following notations are used: S = S+-|-S~ and 
A = S"*" — S~. As we are mainly interested in finding 
the total spin polarization S, A can be eliminated from 
Eqs. ([5]) and ^ via a simple transformation. The final 
equations for three components of electron spin polariza- 
tion in cylindrical coordinates, S = 5*^6,, + Sgeg + 5^6^, 
can be presented in the form 
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where w = w/r is the angular velocity. 

The initial rate of change of spin polarization S may be 
calculated by using the initial condition for A from Eq. 
([5]). In particular, if polarizations of electrons moving 
clockwise and counterclockwise at the initial moment are 
the same, i.e. A(0, t = 0) = 0, then we find from Eq. ([s]) 
that in this case 



as 

~di 



^ 0. 



(10) 



The diffusive limit Eqs. ([7])-(|9| is realized on time scales 
much longer than r. In this case, we can neglect d^Si/dt^ 
terms in these equations and notice that the resulting 
equations are of diffusion type with the diffusion coeffi- 
cient D = v'^T. 



RELAXATION OF HOMOGENEOUS SPIN 
POLARIZATION 

The most interesting and simple solution of Eqs. Q- 
([9| describes the relaxation of the homogeneous spin po- 
larization initially directed along z-axis (see Fig|l]) 



(11) 



where Sq is the initial spin polarization amplitude. Since 
the initial condition (111 is symmetrical with respect to 
rotations about the ring axis, 5^ = at anytime, and 
the spin polarization components Sr and Sz should not 
depend on e. Then Eq. ([8| satisfies identically, and Eqs. 
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(0 and (|9| can be simplified to the following relations 



1 dSr 

T~dt' 
Ids, 



'■Sr + ujns, = 0, 



n^s, + ujnSr = 0. 



(12) 
(13) 



dt2 T dt 

Solving these simple equations with the initial conditions 



( 10 1 and (11), we find that the relaxation of homogeneous 
spin polarization initially directed along the ring's axis is 
described by 
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(14) 
(15) 



where k = y ^ — (w^ + i}'^). It should be emphasized 
that the parameter k can take both real and imaginary 
values depending on system parameters. In particular, 
we can define a diffusive (t^/liP- + fi^ <c 1, Im(K) = 0) 
and ballistic (t\JliP + fi^ ^ \^ Re(K) = 0) limits of spin 
dynamics. According to these definitions, in the diffusive 
limit the mean free path is much smaller than the ring 
radius, £ <C r, and the spin precession angle per mean free 
path is small, r/^ <C 1. In the opposite ballistic limit, at 
least one of the following inequalities should hold: £ ^ r, 
rji^ 1. Any of these two conditions for the ballistic limit 
results in oscillations in spin relaxation dynamics. In 
the diffusive limit, however, the spin polarization decays 
exponentially without any oscillations. 

Eqs. (14 1 and (15 1 can be further simplified in 
the diffusive and ball isti c limi ts. In diffusive limit, 
rVuPm? <^ 1, Eqs. (141 and (15 1 take the form 



Sr{t) 



uiflSo 



1 -e 




(16) 
(17) 



It is clearly seen from the above expressions that in this 
case the relaxation of spin polarization is characterized 
by the time {t{uj'^ + fl^))^^, which is much longer than 
T. In the opposite ballistic limit, when t\/uP~+T}? ^ 1, 



Eqs. (14) and (15) can be reduced to 



Sr{t) 



^2 + fi2 



l-e--^ cos{^/uj^ + nH} , (18) 



re 2t cos(vw 



(19) 



Consequently, the time-dependence of spin polarization 
components involves an exponential decay with a time 
constant 2t modulated by oscillating functions. Fig. [2] 
shows an example of spin relaxation dynamics in both 
diffusive and ballistic limits. 
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FIG. 2: (Color online) Relaxation of electron spin polarization 
in a ring in the (a) diffusive and (b) ballistic regimes. So — 0. 
This plot was obtained using Eqs. ( 14 1 and ( 15 I with the 
following parameter values: fl/ui = —0.5, tk = 0.45 (in (a)) 
and TK — i5 (in (b)). 



The most important feature of the solution given by 
Eqs. (14) and (15) is that the electron spin polariza- 
tion does not decay completely to zero. At long times, 
it shapes into a persistent crown-like spin polarization 
structure (schematically presented in Fig. [ij. The long- 
time values of spin polarization components are 



- 5*0 , 5*2 



;So. 



(20) 



The above relations indicate that during the spin relax- 
ation process the amplitude of spin polarization decreases 
by the factor 



S 
So 



\J S'^ 4" >S'2 
5o 



(21) 



In this persistent crown-like spin polarization structure 
the angle between the ring's axis and the direction of 
spin polarization does not depend on the magnitude of 
the initial polarization Sq and is equal to 



tani/) 



n 

UJ 



-rjr. 



(22) 



It follows from Eqs. (21 1 and (22) that the geometric pa- 



rameters of persistent crown-like spin structure depend 
(besides the ring radius r) only on the ratio of the fre- 
quency of spin-orbit precession to the frequency of the 
electron rotation, Q/uj. Moreover, we note that in the 
case of the opposite direction of the initial spin polariza- 
tion the asymptotic value of Sr changes to the opposite 
one. In particular, if we consider the situation shown 
in Fig. [T| (realizable at a negative value of a), then the 
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change Sq —Sq will result in the opposite directions of 
the radial and z components of spin polarization. Such 
an asymmetry can be considered as a result of the sym- 
metry breaking by the spin-orbit interaction. 




PERSISTENT SPIN STATES 

The crown-like persistent spin polarization structure 
discussed in the previous section is the simplest example 
of persistent spin states of spin kinetic equation (|7])-(|9|. 
In this section we find the complete set of the persistent 
spin states of Eqs. ([7|-([9]). For this purpose, we search 
specific solutions of these equations in the complex form 



Si9,t) 



where n is an integer number. The possible values of 
A and corresponding amplitudes Sr, Sg, Sz have to be 
found. Substituting Eq. ([23| into Eqs. ([7])-([9]) we readily 
obtain a system of linear homogeneous equations. The 
condition of consistency of these equations give rise for 
the following values of A: 



Sr{0,t) \ 






Se{e,t) 


1 o 1 in6 Xt 


(23) 




= ySj j e e , 





A? = 



(24) 
(25) 
(26) 



Solving the linear homogeneous equations for each value 



of A (given by Eqs. (24l-(p6|) we find corresponding 
amplitudes 



^1 




(27) 



It is interesting to note that the vectors given by Eqs. 



(27 1 do not depend on r and n. They are orthogonal and 
normalized (in complex sense) , so it is easy to expand any 
vector (for example, related to initial spin polarization) 
in this basis. 

A spin state is persistent when A = (this follows 
directly from Eq. (23l). It can be seen from Eqs. (24)- 
( 26 1 that there are three possibilities when the condition 
A = is realized. The first possibility is when n — 0, 
in this case, A^ is equal to zero. The corresponding per- 
sistent spin polarization structure is given by S = ASi, 



FIG. 3: (Color online) n — 2 stationary state (as follows from 
Eq. ([29| at 00 = 0). 



that is 



, " A, Sg = 0, , A, (28) 



where A is an arbitrary constant describing the magni- 
tude of spin polarization in this state. It can be seen that 
the persistent spin state given by Eq. ( 28 1 is realized in 



the relaxation of homogeneous spin polarization initially 
directed in z direction (this problem was considered in 
the previous section, see Eq. (20 1). 

As it is seen from Eqs. ( p5| and (26 1, two other pos- 
sibilities (X^ = or Ag" ~ 0) may only be realized if 



02 

1 -h ifj- or n 



' 1 + ^ correspondingly. Since 
n is an integer number, these conditions are fulfilled only 
for certain combinations of system p aram eters as dis- 
cussed below. By substituting Eqs. (27 1 for S2 or S3 
into Eq. (23 1 and taking real or imaginary parts, we ob- 



tain the persistent spin polarization structures that can 
be presented as 



B 

Sr — — sin n 
n 



Sz = 



60), Se ^ Bcosn{e-eo), 
BJv'^ — 1 

— sgn(r!)sinn(0-(?o), (29) 



where n = y 1 + £3- = 2, 3, 4, B and 9q are arbitrary 
real constants, which determine spin polarization ampli- 
tude and angular position of the structure respectively. 
Fig. [3] shows schematically the persistent spin polariza- 
tion structure (29 1 for n = 2. 



It should be noted that the persistent spin polarization 
structures (Eq. ( p9| ) exist only at certain relations be- 
tween the ring's radius r and the strength of spin-orbit 
interaction 77. Since fi/o; = 7]r, the conditions for persis- 
tent spin states are given by 



1, n = 2,3,4, 



(30) 



In addition, we note that at r — >■ cxi, n becomes much 
greater than 1. In this case in accordance with Eq. (30 1 
1, and persistent spin polarization structure 



( 29 1 becomes the usual plain spin helix [F with wave 



vector 77. 
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GREEN FUNCTION 



Linear combinations of specific solutions (|23| (with 
amplitudes (27 1 and corresponding A-s given by Eqs. 
(24l-(26)) of homogeneous linear Eqs. ([7])-(|9]) are also 
solutions of these equations. Therefore, the general solu- 
tion of the system (|7])-(|9| can be presented as a sum 



(31) 



where A'^{n) are complex constants, v — 1,2,3, a — 
+,— , and n runs over all integer values. In order to 
find the constants A1{n) we should specify the initial 
conditions for spin polarization 



[ dt 



S{e,Q). (32) 



Using initial conditions ( 32 1 we obtain 

de^ -rue, (^.,^(go,0)-AT(n)S(go,0)) 







27r 



AJ(n)-Aj(n) 



(33) 



where (</>, ip) — + 4'e''Po + 4'z'4^z is the inner product 
of amplitudes. It is easy to check that {S^, S^) = Su^i- 
The Green function of spin kinetic equation can be 



obtained substituting Eq. (33 I into Eq. The Green 



function can be employed to find spin polarization at any 
moment of time for any given initial conditions as 



dt 



dOoGaii{0 — 9o, t)Si3{6o, 0) 
d9oG^p{0-eo,t)Spieo,O), (34) 



where a or /? take the values r, 0, z and summation over 
repeated indexes is implied. 

In particular, the components of the Green function 
(it is actually 3x3 matrix) of the system ([7|-([9| can be 
written as 



»^Go(0,t)+c.^G,(g,^) 

GrA0,t) - 



G9.(0,t) = 



+ ' 

uj^\G,,[d,t)~G,(6,t)\ 

^2 + 



Gre{0,t) — ^G0r{O,t) 

Gee{0, t) = Gc{0, t), 



(35) 

(36) 

(37) 

(38) 
(39) 
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FIG. 4: (Color online) Spin polarization components at 
— tt/A induced by a delta function excitation at S = 
{Si3ie,t = 0) = 0, 5/3(6', t = 0) = 5{e)6fi,^, where S{..) is the 
delta function and Si^j is the Kronecker delta). This plot was 
obtained using the parameter values rfi = —0.1, tuj — 1 (in 
(a)), and toj = 0.1 (in (b)). 



G..i9,t) = 
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where 
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The above equations can be further simplified. In fact. 



the sums appearing in Eqs. ( 44 1-( 46 1 can be calculated in 
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the following way. First of all, we employ the well-known 
formula from the Bessel function theory 

sin (ty^Lu'^ — q"^^ 
dCe^-«$(t - (g^t^^) , (48) 



\ ^+00 
2 



where ^{t) is the Heaviside step function, and Io{t) is the 
modified Bessel function of zero order. Substituting Eq. 



(48 1 into Eqs. (44l-(46l we perform summation in these 



equations taking into account that 

n— +00 n— +00 



(49) 



n— — 00 



n— — 00 



Then, the integrals in Eqs. (44)-(46l modified by Eq 



(48 1 are trivially evaluated yielding the following results 
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2Tm) 
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In fact, all sums in Eqs. (50)-(52l are finite sums be 



cause of the Heaviside functions $(w< — \9 — 2Tm\). For 
any given values of time t and angle 9 only the terms with 
n in the range [9 — (jjt)/{2iT) < n < {9 + ujt)/{2-K) give 



non-zero contribution to the R.H.S. of Eqs. (50l-(^52 



Fig. |4] shows an example of application of the Green 
function formalism to calculations of spin polarization 
dynamics. Specifically, we plot spin polarization compo- 
nents at = 7r/4 induced by spin polarization excitation 
ai 9 = 0. Such an excitation induces clockwise and coun- 
terclockwise propagating waves of spin polarization that 
finally evolve into a steady homogeneous spin polariza- 
tion. Each time when clockwise and counterclockwise 
propagating waves of spin polarization reach 9 = n/4, 
the spin polarization at this point changes in steps (see 
Fig. |4|a)). The step amplitudes, however, decrease in 



time because of the diffusive component in spin trans- 
port. This type of behavior resembles a multiple echo. 
In addition, in the diffusive limit shown Qb) we did not 
observe clear steps in spin polarization. At long times, 
Se ^ and the crown-like spin polarization configura- 
tion is formed. 



CONCLUSIONS 

Spin relaxation in confined structures exhibits some re- 
markable properties and behavior that are not found in 
infinite 2D systems. In this paper we have studied the 
electron spin relaxation in the ring. It has been found 
that the homogeneous spin polarization along the ring 
axis transforms into a persistent crown-like spin struc- 
ture. Moreover, a family of persistent spin states in the 
ring has been identified. We have also derived the Green 
function of spin kinetic equation and used it to investi- 
gate the propagation of a point spin excitation. 
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